Presnel equations and the refractive index of active media 
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There exists a class of realizable, active media for which the refractive index cannot be defined as 
an analytic function in the upper half-plane of complex frequency. The conventional definition of 
the refractive index based on analyticity is modified such that it is valid for active media in general, 
and associated Fresnel equations are proved. In certain active media, the presence of a "backward" 
wave, for which both phase velocity and Poynting's vector point towards the excitation source, is 
demonstrated. 

PACS numbers: 41.20.Jb, 42.25.Bs 



I. INTRODUCTION 

With the recent advances in the fabrication of pas- 
sive materials and metamaterials P, Q , it seems plausible 
that novel classes of active materials may be developed. 
While the basic electromagnetic concepts, such as micro- 
scopic and relativistic causality, certainly apply to active 
media, the analysis of electromagnetic wave propagation 
may differ. For example, at normal incidence to a bound- 
ary between two passive media, the transmitted field is 
proportional to exp(inwz/c), where n is the refractive in- 
dex, lu is the angular frequency, z is the coordinate in the 
direction away from the boundary, and c is the vacuum 
velocity of light. Since z is arbitrary, relativistic causal- 
ity dictates n to be an analytic function in the upper 
half of the complex w-plane 3, 4]. On the other hand, at 
least in principle there exist realizable, active media for 
which n cannot be identified as an analytic function in 
the upper half-plane, fndeed, consider a Lorentz model 
with relative permittivity e = 1 — /, where / is given by 



icuT 



(1) 



Here, u>q, and T are positive parameters. If F > this 
medium is active; furthermore, if F > 1, e has a simple 
zero in the upper half-plane. Thus, assuming the relative 
permeability /u = 1, when F > 1 the refractive index 
n = contains a branch point in the upper half- 

plane. 

These media would seem to violate relativistic causal- 
ity. This is certainly not the case; to resolve this appar- 
ent paradox, I will provide a rigorous derivation of the 
Fresnel equations in the general case, based on the elec- 
tromagnetic solutions of a finite slab. The definition of 
the refractive index is generalized so that it is valid for 
active media. It will become clear that when the refrac- 
tive index contains branch points in the upper half-plane, 
the Fresnel equations and the refractive index do not 
have unique physical meaning for real frequencies. The 
Kramers-Kronig relations for the refractive index must 
be modified accordingly. 

As an extra bonus, the general framework in this pa- 
per can be used to analyze novel classes of active media. 



Recently, it has been suggested that certain active, non- 
magnetic media with simultaneously passive and active 
material resonances can yield negative refraction I 
will consider the Fresnel equations for these types of ma- 
terials, and demonstrate that at normal incidence, the 
"backward" wave is excited so that both phase veloc- 
ity and Poynting's vector point towards the excitation 
source. While negative refraction is found at oblique in- 
cidence, it is shown that evanescent field amplification, 
similarly to that in the Veselago-Pendry lens [(J 13 > does 
not happen. 

The remaining of the paper is structured as follows: 
In Section II the theory of the electromagnetic parame- 
ters is reviewed. Section III contains an electromagnetic 
analysis of a finite slab surrounded by vacuum. From the 
resulting fields, in Section IV the Fresnel equations are 
proved in the general case, and the associated refractive 
index is identified. The results are interpreted in the form 
of examples (Section V); inverted Lorentzian media, and 
more complex media featuring negative refraction. 



II. THE ELECTROMAGNETIC PARAMETERS 

Any electromagnetic medium must be causal in the mi- 
croscopic sense; the polarization and magnetization can- 
not precede the electric and magnetic fields, respectively. 
This means that the relative permittivity e(w) and the 
relative permeability fi(uj) obey the Kramers-Kronig re- 
lations. In terms of the susceptibilities \ = 6 ~ 1 or 
X = f-i — 1, these relation can be written 



Imx = HRex, 
Rex = -ftlmx, 



(2a) 
(2b) 



where H denotes the Hilbert transform 0, H, 0, Ho| . 
These conditions arc equivalent to the fact that x is an_ 
alytic in the upper half-plane (Imtj > 0), and uniformly 
square integrable in the closed upper half-plane (Titch- 
marsh' theorem) [HHEl. The susceptibilities are de- 
fined for negative frequencies by the symmetry relation 



x(- 



(3) 
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so that their inverse Fourier transforms are real. 

For passive media, in addition to (J2J and © we have: 



Imx(w) > for uj > 0. 



(4) 



The losses, as given by the imaginary parts of the sus- 
ceptibilities, can be vanishingly small; however they are 
always present unless we are considering vacuum ficij| . 
Eqs. ©-© imply that Imx > in the first quadrant 
(Reui > and Imw > 0) , and that 1 + x 1S zero- free in 
the upper half-plane [lfjj. Thus for passive media, the 
refractive index n = ^fe^/Ji can always be chosen as an 
analytic function in the upper half-plane. With the ad- 
ditional choice that n —> +1 as w -> oo, rt is determined 
uniquely by 



VRR exp [i (arg e + arg fi)/2] 



(5) 



where the complex arguments are restricted to the inter- 
val (— 7r,7r]. From Eq. (JSJ it is immediately found that 
© and (@J hold for the substitution x —> n - Moreover, 
the Kramers-Kronig relations for the refractive index are 
established using Titchmarsh' theorem by noting that 
n — 1 is analytic in the upper half-plane, and that n — 1 
has the required asymptotic behavior: The fact that x 
satisfies means that x satisfies the square integrabil- 
ity condition. From n = y/e^/JI, it follows that n — 1 has 
a similar asymptotic behavior. This gives 



Ivan — TC {Ren — 1}, 
Re n — 1 = — TL Im n. 



(6a) 
(6b) 



On the basis of causality, it is clear that the suscepti- 
bilities of active media also satisfy J2J an d ©■ However, 
since Q does not hold for active media, it is no longer 
true that e and \x are zero-free in the upper half-plane. 
Indeed, we have already seen that the permittivity asso- 
ciated with an inverted Lorentz model may have a simple 
zero in the upper half-plane. Thus, the refractive index 
cannot always be chosen as an analytic function in the 
upper half-plane. 

To clarify this point, I will consider one-dimensional 
electromagnetic wave propagation in an active medium. 
The main goal is to identify the reflection and transmis- 
sion at a single interface, and to prove that these waves 
obey relativistic causality. The Fresnel equations could 
be proved in the conventional way; however, it might not 
be clear a priori how to specify the wavevector of the 
transmitted wave in the active medium. In fact, con- 
trary to the case with passive media, it is not necessarily 
correct to let the Poynting vector of the transmitted wave 
point away from the interface. Moreover, it is not always 
possible to choose the wavevector as an analytic function 
for Imcj > 0. Therefore, the active medium is assumed 
to have finite thickness initially. Then the transmitted 
wavevector at the far end of the slab can be specified 
trivially, and the electromagnetic solutions can be found. 
Subsequently, using the result of the finite slab, the fields 
of a half-infinite medium can be determined. 



slab 



z=0 



z=d 



FIG. 1: A plane wave incident to a slab of thickness d. 



III. LAPLACE TRANSFORM ANALYSIS OF AN 
ACTIVE SLAB 



Consider a plane wave in vacuum, normally incident 
to a linear, isotropic, and homogenous slab with param- 
eters e and fi, see Fig. ^ While vacuum is chosen as 
the surrounding medium, it should be clear that any 
other passive media yield similar results. The slab is 
located in the interval < z < d, where d is the thick- 
ness of the slab. A causal excitation on the left-hand 
side of the slab (z = 0~) is assumed, i.e., the fields for 
z > vanish for time t < 0. Since the medium is active, 
the physical, time-domain fields may diverge as t — > oo, 
so Fourier transformed fields do not necessarily exist. 
(Strictly speaking, in practical materials, the fields do 
not diverge. When the fields become sufficiently strong, 
they deplete the gain. However, we restrict ourselves to 
the linear regime of the susceptibilities, and bear in mind 
that the model breaks down when the fields become large. 
The divergences are indicators of instabilities such as, for 
example, the ignition of lasing.) Thus, instead of Fourier 
transforms, we use Laplace transforms of the fields. De- 
noting the physical, time-domain electric field £(z, t), the 
Laplace transformed field for z > is 



/>oo 

E(z,uj) = / £(z,i)exp(i 



(7) 



where Imw > 7, for a sufficiently large, real constant 7. 
Assuming that the source is located at z = z s < 0, Q 
still applies for z s < z < provided the lower limit in 
the integral is replaced by z s /c. Since the fields every- 
where are assumed to vanish for t equal to the lower limit 
in the integral, the partial time derivatives in Maxwell's 
equations become simply — iu in the transform domain. 
At the boundaries the Laplace transformed electric and 
magnetic fields must be continuous. Thus we can use a 
similar analysis method as in the Fourier case, consid- 
ering the transformed fields at each to. Assuming the 
forward propagating electric field £ + (z,t) for z < 0, and 
recalling that £(z,t) for z > d is assumed to vanish for 



3 



negative time, we can write 



IV. FRESNEL EQUATIONS 



E(z,uj)/E+(0,co) (8) 

exp(iuiz/c) + Rexp(—iuz/c) for z s < z < 0, 
= I S + cxp(ikz) + S~ exp(-ikz) for < z < d, 



T exp(iujz/c) 



for z > d, 



where E + (0, u>) is the Laplace transform of the excitation 
£ + (0,t). By matching the electric and magnetic fields at 
both interfaces, the total reflection coefficient R, field 
amplitudes in the slab (S + and S~), and transmission 
coefficient T are found to be 



R = 



T = 



(rf — 1) cxp(— ikd) — (rj 2 — 1) exp(ikd) 
(77 + l) 2 cxp(— ikd) — (77 — l) 2 exp(ikd) ' 

277(77 + 1) 

" (r? + l) 2 - (t? - l) 2 exp(2ikd) ' 

277(77-1) 

" (r? - f ) 2 - (77 + l) 2 exp(-2ikd) ' 
477 

(7/ + f ) 2 exp(— ikd) — (77 — l) 2 exp(ifcc?) 



(9a) 
(9b) 

(9c) 
(9d) 



Here, k = noj/c, 77 = [ijn, and n — y/ejl. Note that R, 
S = S + cxp(ikz) + S~ exp(-ikz), and T are unchanged 
if 71 — * —77, so the choice of the sign of 71 does not matter. 
In other words, only even powers of n are present in R, 
S, and T, which means that the fields contain no branch 
points in the upper half-plane. 

For active media, R, S, and T may contain poles in the 
upper half-plane. Despite poles, the fields arc relativisti- 
cally causal. This fact is established using the analyticity 
of the respective functions in the half-plane Im u> > 7 and 
their asymptotic behavior in that region. For example, 
T tends to exp(iuid/ c) as Imw — > 00, which means that 
the signal- front delay through the slab is d/c [T3 |. 

The time-domain solutions are found by inverse 
Laplace transforms of the associated fields. For example, 
if the excitation is £ + (0, t) = u(t) exp(— iuiit), where u(t) 
is the unit step function and Imwi < 0, the time-domain 
field in the slab is 113 



£(z,t) 



1 

2^ 



27 + 00 



S 



iu>i — ioj 



■ exp(—iujt)dtjj. (10) 



When S has no poles in the closed upper half-plane, 7 
can be set to zero. It follows then that the field S can be 
interpreted in the usual way at real frequencies. When S 
contains poles in the upper half-plane, the integration 
path in (|10() should be located above all of them. If 
we still insist on locating the integration path along the 
real axis, the resulting discrepancy would correspond to 
the sum of residues of the poles in the upper half-plane. 
Such residues are exponentially increasing and demon- 
strate clearly the fact that when S contains poles in the 
upper half-plane, the slab is electromagnetically unsta- 
ble. 



To facilitate the analysis and interpretation of active 
media, it is useful to obtain the fields when the slab fills 
the entire half-plane z > 0. For passive media, one can 
take the limit d — > 00 in JjJ) to obtain 



R = 



?7-l 
77 + I' 



S = 



2rj_ 
77 + I 



exp(inujz/c). 



(11a) 
(lib) 



In obtaining Ijllfl the sign of n has been chosen such 
that Imn is positive in the first quadrant (Re a; > and 
Im uj > 0). For passive media, this choice is equivalent 
to the choice in (J5J , or in the conventional definition of 71 
based on analyticity. Eqs. I jlip are essentially the Fresnel 
equations; however the propagation factor exp(i7io;2:/c) 
has been inluded in the transmitted field, to emphasize 
the position dependence of the field in the active medium. 

For active media, one could guess that (jllfl remains 
valid, but it is not clear a priori how to choose n. The 
imaginary parts of e and [i may take both signs in the 
first quadrant, and it may not be possible to choose n 
as an analytic function in the upper half-plane. Sur- 
prisingly, as shown by the following example, taking the 
limit d 00 in Q leads in general to unphysical re- 
sults. Consider a slab of an inverted Lorentzian medium 
e = 1 — / and /i — 1, where / is given by l[T]l. Assuming 
F is small (corresponding to low gain), n w ±(1 — J/2). 
In the limit d — > 00, we find R = (77 + l)/(?7 — 1) and 
S = 2?7exp(— inu)zjc]j(r\ — 1), where n = 1 — //2. This 
solution is clearly unphysical as i?, 5 — > 00 as w -> 00; a 
convergence half-plane Imw > 7 does not exist. 

To find the correct solution, we start with the time- 
domain solution i|10[) . where d is finite. Due to the finite 
bandwidth of physical media, |x| cx 1/ui 2 as ui — > 00 in 
the upper half-plane Thus, for a sufficiently large 
\i>}\, e w 71 w 1. Choosing the sign of n such that 
n ~ 1 here leads to \n — 1| cx 1/lo 2 for uj — > 00. 
This implies in turn that there exists a 7 > such 
that I [(77 — 1)/ (77 + l)]exp(ikd)\ < 1 for any oj with 
Imw > 7. We can now expand (|9bl) and l(9"cl) into ge- 
ometrical series in [(77 - l)/(?7 + l)] 2 exp(27:fcd). Hence, 
S = S + exp(ikz) + S~ exp(-ikz) becomes 



S = 



2rje tkz 
77 + I 



77-I 



. 2ik-cl 



(12) 



277(77 — l)e 



likd—ikz 



(v + 1) 2 



1 



Jlikd 



After substitution of S into 1|1(J[I . and using the signal- 
front delay property of the Laplace transform [14], it is 
apparent that only the first term leads to a nonzero re- 
sult for t < 2d — z. This term coincide with (jllb|) . In a 
similar way, we find that for time t < 2d, the reflected 
wave at z — 0~ is given by the term (jllafl substituted 
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for S in i|10|) . Now we can clearly take the limit d — > oo. 
This will not alter Qllf). but will ensure the validity of 
the associated time-domain solutions for all times. The 
solution has now been found; however, the evaluation of 
the fields far away in the upper half-plane is imprac- 
tical. Fortunately, by Cauchy's integral theorem we can 
move the line 7 down towards the first branch point of n 
or zero of r\ + 1 . If there are no such points in the closed 
upper half-plane, we can even evaluate on the real 
axis. Note that although the original functions R and 
S from Q are not analytic in the upper half-plane, the 
terms (1111) are analytic provided r; + 1 has no zeros and 
n is analytic. Also note that in all cases, the expressions 
(fTT|) are analytic in some region Imw > 7, and S tends 
to exp(iu;z/c) as Imw — > 00. In other words, £ (z,t) = 
for t < z/c [14], and relativistic causality is guaranteed. 

We conclude that the Fresnel equations l|ll|) retain 
their form, but not always their interpretation for active 
media. When e/i has no odd-order zeros in the upper half- 
plane, the sign of the refractive index is determined such 
that n^+1 as ij ^ 00, and such that n is analytic in 
the upper half-plane. Furthermore, if 77 7^ — 1 everywhere 
in the closed upper half-plane, the Fresnel equations l|ll|) 
can be evaluated and interpreted in the usual way along 
the real frequency axis. On the other hand, if e(i has odd- 
order zeros, or r\ + 1 has zeros in the upper half-plane, 
(|llfl contain points of non-analyticity. The reflected and 
transmitted fields (|llf) should then be evaluated along 
the line Imw = 7, where 7 is larger than the maximum 
Im lo of the points of non-analyticity. The sign of the re- 
fractive index is determined so that n — > +1 as w - > 00, 
and so that n is analytic for Imu > 7. Note that when 
efj, has odd-order zeros in the upper half-plane, n is not 
uniquely determined along the real frequency axis, and 
consequently, one must be careful with its interpretation 
there (See Appendix A). 

For non-orthogonal incidence we can use a similar ar- 
gument as that given above. When the incident wave 
vector is k = (k x ,k y ,k z ), the TE fields of a finite slab 
are given by the substitutions 77 — > (ik z /k' z and k — > k' z 
in ©. Here, k' z = e(iuj 2 /c 2 — k 2 — k 2 . Again, we note 
that the choice of sign of k' z does not matter as only even 
powers of k' z are present. Following the argument leading 
to (|llfl for active media, we arrive at the similar Fresnel 
expressions 



S=-^_exp(< 2 ), (13b) 

valid for Imw > 7. Here, the sign of k' z is determined 
such that k' z is analytic for Imw > 7, and such that 
k z — > +lu/c asw-^oo. The parameter 7 is chosen larger 
than the maximum imaginary part of branch points of k' z 
or zeros of (ik z + k' z . If there are no such points in the 
closed upper half-plane, 7 can be chosen to be zero, i.e., 
the Fresnel equations can be interpreted along the real 
frequency axis in the usual way. 



Since the refractive index is analytic in the half-plane 
Imu > 7, and since n — 1 has the required asymptotic 
behavior for Imu > 7, the Kramers-Kronig relations for 
active media can still be written as in JSJ; however, the 
integrals are taken along the line Imw = 7. Note that 
the conventional Kramers-Kronig relations for n, as in- 
tegrated and evaluated along the real frequency axis, are 
not valid. 



V. EXAMPLES 

A. Inverted Lorentzian media 

In this example, we consider an inverted Lorentzian 
medium with e = 1 — /, where / is given by (T). and 
(1 = 1. Since / ^ everywhere, rj + 1 has no zeros. Con- 
ventional optical gain media, such as the Erbium-doped 
silica fiber, have relatively low gain per length of fiber. 
We will therefore first consider the well-known case where 
/ satisfies |/| <C 1. Then e has no zeros, and following 
the recipe above, we find that the Fresnel equations ljll|l 
can be interpreted straightforwardly for real frequencies, 
and n « 1 — //2. 

Lifting the assumption of low gain, we observe 
that the zeros of e are located at uj± = — iT/2 ± 
i v /r 2 /4+ (F- l)u 2 . When F > 1, the zero u+ is lo- 
cated in the upper half-plane. Consequently, in this case 
neither n nor the Fresnel equations have unique physical 
meaning along the real frequency axis. The fields should 
instead be evaluated along the line Imw = 7, where 7 is 
chosen such that 7 > Im lo + . 

It is useful to recall that the fields, as given for 
Imu > ImW-)., can be interpreted by choosing an ex- 
ponentially increasing excitation. For example, we may 
choose an excitation £ + (Q.t) — u{t) exp(— iu)\t), where 
Imu + < Imui < 7. For sufficiently large t, the time- 
domain field HI Up gets the main contribution from the 
residue of Sexp(— iu>t)/(iu>i — iu>) at w\. In other words, 
S is roughly the complex amplitude of the exponen- 
tially increasing wave. The resulting field dependence 
exp(inu>iz/ c — iuiit) means that the wave has its "phase 
velocity" in the positive z-direction. (For this medium 
Ren > in the first quadrant of complex frequency.) 

The exponentially increasing excitation is somewhat 
artifical, and it would be interesting to evaluate the time- 
domain solution when Imct>i < 0, i.e., when the excita- 
tion pulse envelope is exponentially decreasing. In Ap- 
pendix A, the time-domain field at a fixed z is evaluated, 
and it is argued that the pulse diverges as t — » 00. Note 
that the field in the medium diverges as a result of the 
medium alone, not as a result of, for instance, amplified, 
multiple reflections. 
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B. Excitation of a "backward" wave 

Consider the medium e = (1 + f) 2 and fi= 1, where / 
is still given by (I). The associated susceptibilities satisfy 
P)l and J3J|, so at least in principle, this medium is realiz- 
able. The product e/i is zero-free in the upper half-plane, 
and we can determine the refractive index such that it is 
analytic in the upper half-plane (and such that n — > +1 
as uj — ► oo). This gives n = 1 + / and r/ = 1/(1 + /). 
Note that 77 7^ — 1 for Imw > 0. It follows that the Fres- 
nel equations 1)11(1 can be interpreted straightforwardly 
along the real frequency axis. 

When F is sufficiently large, there exists a real fre- 
quency u>i = u)oy/l + F/2 with Ren(wi) « —1 while 
Im?i(wi) w ■^-^/l + F/2 is small. This seems to be a 
potential route to negative refractive index at optical fre- 
quencies Since Ree(cji) ~ 1 and /x = 1, not only the 
phase velocity but also the Poynting vector point in the 
negative z-direction at this frequency. As lm.n(u>i) > 
and Ime(wi) < 0, the field can be interpreted as a grow- 
ing wave propagating in negative z-direction. Neverthe- 
less, the excitation of this "backward" wave does neither 
rely on finitencss of the thickness of the active medium 
nor excitation from z = +00. Note that any causal exci- 
tation involves an infinite band of frequencies, and out- 
side the limited negative-index band, the phase velocity 
and Poynting vector point in the forward direction. 

For this medium, the transmitted time-domain field 
at z = + can be calculated analytically. In this 
case we choose a real excitation of the form £ + (0,t) = 
u(t) cos(wii), where u>i = wo\A + F/2. With the help 
of the corresponding Laplace transform and the Fresnel 
transmission coefficient 2/(n+l), we find after some sim- 
ple algebra: 



15 



27+00 



(bj 2 + iu)T — lOq)lu e~xjp(—iu>t) 



du>. 



£ ( 0,t } 2 7 ry i7 _ 00 (u - w 2 )(w - uj-2){u 2 

(14) 

Here 7 = + , and lj±2 — ±tc>i — iT/2 to first order in 
T/u>i. The inverse transform (|14f> can be found exactly 
by calculating the residues in the closed, lower half-plane. 
Assuming T/u>i <C 1, this procedure yields 



£(o,t) 



cos(uit), 



Im n{ui 



for < Tt/2 <C 1, 
■sin(wit), forrt/2>l. 



(15) 



(The steady-state solution could certainly be found di- 
rectly from l|llb|l .l We observe that the transmitted 
field is equal to the excitation for small t. Thus, initially 
the reflection from the boundary vanishes. On the other 
hand, when t is large, the field has grown considerably. 
This fact is interpreted by noting that the "backward 
wave" , with phase velocity and Poynting's vector in the 
negative z-direction, has been built up. This wave draws 
energy from the active medium, yielding a large field at 
the boundary. Note that also the "reflected wave" on the 
left-hand side of the boundary now is large. Nevertheless, 
the fields are produced in a stable manner. A plot of the 
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FIG. 2: The electric field £(z,t) of the example in Subsection 
IV Bl The frequency loq = 1 (normalized), and t = 30. The 
parameters used are F/luo = 0.1 and F — 10. The field was 
calculated by a straightforward numerical evaluation of 11411 
using 7 = O.OIojo- The field plot for small z shows the "back- 
ward wave" , with phase velocity and Poynting's vector in the 
— z-direction. The forerunner, located close to z = 30c, moves 
at a speed c in the +z-direction. 



backward wave and the associated, forward propagating 
forerunner is given in Fig. [21 

At oblique incidence the backward wave is refracted 
negatively. Although this may seem obvious intuitively, 
the fact can be proved by calculating k' z and identifying 
its sign along the real frequency axis using the proce- 
dure outlined above. Despite negative refraction, the 
medium has rather different properties from the left- 
handed negative-index medium considered by Veselago 
and Pendry Considering a finite slab of sufficiently 

small thickness such that T has no poles in the upper half- 
plane, we find T w exp(+zo;id/c) at the frequency u>%. 
Also, evanescent-field amplification, as in the Veselago- 
Pendry lens, will not arise in this medium. For TE polar- 
ization this is seen by the substitutions 77 — > fj,k z /k' z and 
k -> k' z in pd). As before, k' 2 = efiuj 2 /c 2 -k 2 - k 2 . For 
an incident evanescent field with sufficiently large k y , k' z 



is close to an imaginary number. Substitution into (|9c 
gives exponential decay. Note that the choice of sign of 
k' z does not matter. 

An interesting point arises by comparing the results 
in this example with those of the previous example. As- 
suming low gain in the previous example, e w 1 — ia for 
some real frequency uj\. Here, a is a small positive num- 
ber. In the present example, at the frequency u>i we can 
also write e ~ 1 — ia for a small positive a. Indeed, it is 
possible to tune the parameters in the two examples such 
that the two e's (and fi's) are identical at this frequency. 
Nevertheless, the Fresnel equations (|TT)l give completely 
different answers; in the previous example the "forward" 
wave, with phase velocity and energy growth in the -In- 
direction, is excited, whereas in the present example a 
"backward" wave is excited. The dilemma is resolved 
by noting that although the e's are identical for a single 
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frequency, they are not identical globally. As a result, 
the causal excitation at z — — , which necessarily con- 
tains an infinite band of frequencies, will have different 
consequences in the two situations. 



C. Instability at the boundary 

In the previous example, although n ss — 1 at the real 
frequency cji, n = — 1 is exactly fulfilled only in the lower 
half-plane (at the complex frequencies In principle, 

one can construct a causal, active medium for which n = 
— 1 at a real frequency or even in the upper half-plane. 
Indeed, let e = (1 + g) 2 and \i = 1, where 



iGlo 



iloT 



(16) 



We assume T/u -C 1 and G ~ T. The associated suscep- 
tibility x = e — 1 is clearly analytic and uniformly square 
integrable in the closed upper half-plane; thus it satis- 
fies the Kramers-Kronig relations jlfij . By determining 
the refractive index such that it is analytic in the upper 
half-plane, we obtain n = 1 + g. Thus we can find the 
frequencies u> a ± where n = a is fulfilled: 



2u. 



a± 



±2w -iT-iG/(a- 1). 



(17) 



Here we have assumed that \a— 1| > 1 and omitted second 
or higher order terms in T/lu . Using l|17|) we find that if 
G > 2T, the solutions to the equation n = — 1 are located 
in the closed, upper half-plane (equality implies locations 
on the real axis). 

Substituting Ijllbjl into (jl()|l , it is realized that the lo- 
cations of the zeros of n + 1 are crucial for the electro- 
magnetic stability. Clearly if l/(n + 1) has poles in the 
upper half-plane, the associated residues when calculat- 
ing l|l(J|) grow exponentially with t. In other words, the 
system consisting of the two media (vacuum and the ac- 
tive medium) is electromagnetically unstable. For a given 
active medium, this instability may be eliminated if the 
vacuum on the left-hand side of the boundary is replaced 
by another passive medium: Let the new medium have 
fi = 1 and (approximately) constant refractive index n p 
for | lu — ujq | < r. The instability would not be present 
provided n ^ — n p in the upper half-plane. According to 
(fT7|l . if 2T < G < (rip + l)r, there is no instability, in 
contrast to what we found above for vacuum. 



VI. CONCLUSIONS 

General forms of the Fresnel equations, valid for ac- 
tive media, have been proved. When the refractive index 
contains branch points in the upper half-plane, the in- 
dex is not determined uniquely along the real frequency 
axis, but is dependent on the choice of branch cuts. Thus 
the Fresnel equations and the refractive index cannot be 



interpreted along the real frequency axis unless their be- 
havior along the branch cut are taken into account as 
well. The expressions should rather be evaluated or in- 
terpreted in some upper half-plane ImcJ > 7, where the 
line Imw = 7 is located above the non-analytic points. 
The presence of branch cuts in the upper half-plane leads 
generally to instabilities in the sense that any vanishing 
small pulse excites an infinite field (or saturation) in the 
medium. 

The sign of the refractive index should be defined using 
the asymptotic form of n (n — > +1 as lu — > 00), and ana- 
lyticity for Imw > 7. It would have been useful to have 
an explicit relation for the refractive index, similar to JSJ. 
However, for active media such a relation is not valid as 
the complex arguments of e and \x may take any values. 
Nevertheless, if 7 is above all non-analytic points of n, n 
is analytic for Imw = 7. Thus the sign can be identified 
up to a global sign by unwrapping the phase arg e + arg fi 
before using J5j. The global sign is determined easily by 
requiring n — > +1 as to — > 00. 

For most media of interest, there are no branch points 
of 11 for Imcj > 0, and the sign of n is to be determined 
along the real frequency axis. In principle, n is not neces- 
sarily analytic for real frequencies (since e and fi are not 
necessarily analytic there). With the extra assumption 
that efj, is continuous for real frequencies, e/z is continuous 
for Imw > 0. Then the sign of n may still be found by 
unwrapping the phase arg e + arg /j, in (JSJ) , and requiring 
n — > +1 as lu — > 00. 

Using the general framework, we have seen that in cer- 
tain active media, only the "backward" propagating wave 
is excited. For this wave both phase velocity and Poynt- 
ing's vector point towards the excitation source. These 
media are conveniently described using a negative refrac- 
tive index. Nevertheless, evanescent-field amplification, 
similar to that in the Veselago-Pendry lens, does not ex- 
ist. 



APPENDIX A: TRANSMITTED FIELD WHEN n 
HAS BRANCH POINTS 

I will here discuss the time-domain field i|10|) associated 
with the solution (|llbl) . in order to provide an interpre- 
tation of media where n has branch points in the upper 
half-plane. For simplicity the medium is chosen to be of 
the type discussed in Section TV Al The electromagnetic 
parameters are e = 1 — / and /x = 1, where / is given by 
(JTJ, and F > 1. This leads to a single branch point of n 
at lo = uj + , where Rcw + = and Imw + > 0. 

The branch cut is chosen along the imaginary fre- 
quency axis, from the branch point at lo + towards — ioo, 
see Fig. [3] With this choice, the time-domain field can 
clearly be calculated by locating the Bromwich path in 
(fTUfl above lo + (7 > Imw+): 

1 r l i+°° exp(inujz/c) exp(— iuit) 

£(z,t) = — j ; a w . ; ; QUI. 
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FIG. 3: Branch cut and integration paths. The Bromwich 
path Imoj = 7 can be deformed into a path consisting of two 
contributions: Integration along the real frequency axis and 
integration around the branch cut. 



It would however be interesting to locate the integration 
path along the real axis, to investigate the interpretation 
of the refractive index for real frequencies. Then, to ob- 
tain a correct result, we must add an integral around the 
branch cut. That is, 



£(z,t) =£j?(z,t) + A£(z,t), 



(A2) 



where £p{z,t) is the contribution from the integration 
along the real axis (inverse Fourier transform), and 



A£(z,t) results from the integration around the branch 
cut (see Fig|3J). Since n is bounded for Imu > 0, the con- 
tribution from the integration along the circle centered 
about a;+ can be neglected when the radius approaches 
zero. Thus we obtain 



A£(z,t) 



exp(w^) 



(A3) 



cxp(ncjiz/c) exp(— nujiz/c) 



1 



1 



duj. 



where n = n{iuji) is evaluated along the right-hand side 
of the branch cut. For example, considering small z such 
that \nujiz/c\ <C 1 along the branch cut, it is clearly seen 
that A£{z, t) — > oo as t — > oo. Since the integrand in 
I|A1(I is square integrable along the real axis, so is £p(z, t). 
Consequently, the electric field £ (z, t) diverges as t — ► oo. 

The branch cut in Fig. |3] could of course be chosen 
in a different manner, as long as it is located under the 
line Imw = 7. For each choice of branch cut, the sum 
£p(z, t) + A£(z, t) will be the same; however the two sep- 
arate terms will generally be different if the branch cut 
crosses the real axis at a different point. This shows that 
the inverse Fourier transform part £?{z,t) has no phys- 
ical meaning unless it is seen together with A£{z,t). In 
other words, the refractive index has no unique, physi- 
cal interpretation along the real frequency axis unless its 
properties along the branch cut are taken into account. 
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